Las formulaciones integrales son metodologías que consisten en transformar los operadores diferenciales a operadores integrales incorporando de manera exacta las condiciones de frontera. Este trabajo se extiende el uso de formulaciones de ecuaciones integrales (FEI) para el desarrollo de esquemas numéricos para sistemas generales tipo reacción-difusión. Debido a la propiedad de incorporación exacta de las condiciones de frontera, los esquemas numéricos resultantes no requieren de aproximaciones en las fronteras, lo que conduce a esquemas numéricos con ordenes globales de aproximación de O(h²). Las FEI se caracterizan por su metodología sistemática y su estructura matemática que permite la fácil interpretación física de los procesos involucrados en los sistemas reacción-difusión. Para evaluar la capacidad de aproximación numérica del esquema propuesto se consideraron diferentes condiciones que surgen en el estudio de las pastillas catalíticas. Las simulaciones numéricas realizadas muestran que las FEI exhiben mejores aproximaciones numéricas que los esquemas clásicos de diferencias finitas (DF).